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Abstract. In this article we study the dynamics generated by 
germs of parabolic diffeomorphisms / : (C,0) — » (C, 0) tangent to 
the identity. We show how formal classification of a given parabolic 
diffcomorphism can be deduced from the asymptotic development 
of what we call directed area of the e-neighborhood of any orbit 
near the origin. Relevant coefficients and constants in the develop- 
ment have a geometric meaning. They present fractal properties of 
the orbit, namely its box dimension, Minkowski content and what 
we call residual content. 
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1.1. Motivation. Each germ of a parabolic diffeomorphism in the 
complex plane, 



(1) f(z) = z + ai z k+1 + a 2 z k+2 + ..., keN, di G C, ai ^ 0, 



can, by formal changes of variables, be reduced to the standard formal 
normal form 



for an appropriate choice of k G N and a G C. The formal type of a 
parabolic diffeomorphism is given by the pair (k,a), keN, a G C. 
The article is motivated by the following problem: 

Problem 1. Can we recognize a diffeomorphism by looking at one of 
its orbits? 

More precisely, the idea is to read the formal type of a diffeomor- 
phism from the fractal properties of any orbit. 

By fractal properties, we mean box dimension and Minkowski con- 
tent of the orbit, which by definition are computed from the rate of 
growth of the area of ^-neighborhoods of the orbit, as e — > 0. More 
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generally, we can refer to the asymptotic development of this 
a fractal property of the orbit. Using only the area, we noticed that 
the answer to the above question is negative. From the asymptotic de- 
velopment of the area of an e-neighborhood, only information on real 
part of a G C can be obtained. Therefore we first generalize the notion 
of area of the ^-neighborhood of the orbit to be a complex number 
whose modulus is the area and whose argument refers to the direction 
of the orbit in the plane. We call it the directed area. We show in 
the article that formal type can be read from two coefficients and the 
leading exponent in the asymptotic development of the directed area 
of the ^-neighborhood of any orbit, as e — > 0. 

In some sense this problem is similar to the famous problem about 
hearing the shape of a drum, see Section |5j 

Let us comment on applications. We show that the directed area of 
the e-neighborhood of any orbit has an asymptotic development in the 

1 I 1 1—1— — 1 J k n 

scale {e , e , . . . , e , £ loge, . . .}. One can study numeri- 
cally the directed areas of ^-neighborhoods of just one orbit, for small e. 
By comparing them to the scale above, one obtains relevant coefficients 
and concludes the formal normal form of the diffeomorphism. 
It is natural to ask now the converse question. 

Problem 2. If we only know the formal type of a given diffeomorphism, 
can we uniquely determine box dimension and Minkowski content of its 
orbits ? 

Equivalently, we can ask if all the diffeomorphisms inside one for- 
mal class have the same fractal properties. It turns out that the box 
dimension is invariant by the changes of variables inside the formal 
class. This is not the case for Minkowski content and what we call 
residual content if we work with general formal changes of variables. 
However, the problem is solved if we restrict the definition of formal 
equivalence relation and allow only formal changes of variables that 
are tangent to the identity. In this sense, each parabolic germ Q is 
formally conjugate to a simpler germ 

(3) g (z) = z + a lZ k+1 + a\-a- z 2k+1 , 

where a\ G C is the first coefficient of the initial diffeomorphism ([!]) 
and a, k are as in (|2]). We call ^ the extended formal normal form. 
The formal type of a diffeomorphism is thus not completely described 
by the pair (k, a), but by the triple (k, ai,a). 

Finally, with this restricted notion of formal equivalence, we get a 
bijective correspondence between the formal type of a diffeomorphism 
and the fractal properties of its orbits. We comment on the prospects 
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and problems of analytic classification of diffeomorphisms using fractal 
properties of their orbits in Section [5j 

In the real similar idea that fractal properties of an orbit 

near a fixed point, i.e. the rate of growth of its e-neighborhoods, carry 
some information on the properties of the generating function itself, 
was discussed before in e.g. [9] and [H]. A bijective correspondence 
was found between the multiplicity of a fixed point of a function on the 
real line and the rate of growth of ^-neighborhoods of any orbit. 



1.2. Definitions and notations. Let us recall precisely the main def- 
initions and notations we use in this article. 

Let / : C — > C, / G Diff(C,0), be a germ of a diffeomorphism 
fixing the origin. We say that the germ / is parabolic if /'(0) = 1. If 
/'(0) = exp(27riA), where A G Q, the diffeomorphism can be reduced 
to the previous case by considering its higher iterates, but we will not 
discuss it in this article. Therefore, in the neighborhood of the origin, 
we suppose that f(z) is of the form 

f{z) = z + a lZ k+1 + a 2 z k+2 + a 3 z k+3 + o(z k+3 ), 

where a\ G C, a, G C, % = 2, 3, . . . and k G N, k > 1. 

By Sf(zo), we denote the orbit generated by f(z) with the initial 
point z in the neighborhood of the origin, S^(zq) = {z n = f on (z )\n G 
No}. Near the origin, such orbits form the so-called Leau-Fatou flower, 
see e.g. [TO] or [8]. In short, there exist k attracting and k repelling 
sectors, called petals, around equidistant repelling and attracting di- 
rections. Attracting and repelling directions are normalized complex 
numbers (— ai) -1 ^, a l 1 ^ k respectively. Orbits are tangent to attracting 
or repelling directions at the origin, see Figure L2 In the sequel, we 
suppose that zq belongs to an attracting sector of the origin. Otherwise, 
if Zo belongs to a repelling sector, we consider the inverse diffeomor- 
phism / _1 (z) instead. 

Now we discuss some properties of measurable sets in the plane. Let 
U C I 2 , or [/ C C, be a measurable set whose center of mass is not 
the origin. By A(U), we denote its area. In Definition [I] in Section [2j 
we define the directed area A (U) of U as the complex number which 
encodes the area of the set U, as well as the direction in which the set is 
placed. Note that the directed area does not verify the finite stability 
property, that is, A C (U U V) = A C (U) + A C (V), for disjoint sets U 
and V. Thus it does not satisfy the properties of a vector measure 
defined in e.g. [5]. Moreover, this notion should not be confused with 
directional e-neighborhood, also called directional Minkowski sausage, 
and directional Minkowski content, defined in 
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FIGURE 1. Attracting and repelling sectors and directions for e.j 
f(z) = z + z 4 + o(z i ). 



The fractal properties of a set U are related to the asymptotic behav- 
ior of the area of its e- neighborhood, denoted U £ , as £ — > 0. The growth 
rate of the area A(U £ ) of the ^-neighborhoods reveals the density of the 
set. It is measured by box dimension and Minkowski content of U. We 
recall the definitions of the Minkowski content and the box dimension 
of a measurable set [/cK 2 (or C). 

By lower and upper s- dimensional Minkowski content of U, < s < 
2, we mean 



Ml{U) 



lim inf 

£->0 



A(U £ 



and M* S (U) = lim sup 



A(U £ 



£-5-0 



respectively, 
defined by 



Furthermore, lower and upper box dimension of U are 



diYn P U = inf{s > | M%{U) = 0}, dim B U = mf{s > | M* S {U) = 0}. 

As functions of s e [0,2], M* S {U) and Ml(U) are step functions that 
jump only once from +oo to zero as s grows, and upper or lower box 
dimension are equal to the value of s when jump in upper or lower 
content appears. 



If dim p U 



dim^t/, then we put dim^(t/) 



dim B U 



dini£?7 and 



call it the box dimension of U. In literature, the upper box dimension 
of U is also referred to as the limit capacity of U, see [TT] . 

If d = dim B (U) and, moreover, Mf(U) = M* d {U) G (0, oo), we say 
that the set U is Minkowski measurable. In that case, we denote the 
common value of the Minkowski contents simply by Ai(U), and call it 
the Minkowski content of U . 
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In short, if A(U e ) ~ Ce 2 s , as e -> 0, for some s e [0,2], C > 0, in 
the sense that lim^o ^2-^ = C> then dim^(?7) = s and C/ is Minkowski 
measurable with Minkowski content M.(U) = C. For more details on 
box dimension, see Falconer [3] or Tricot j!3j . 

In this article, previous definitions are considered for an orbit of a 
parabolic diffeomorphism which accumulates at the origin, U = S*(zq). 
The asymptotic behavior of the directed area of its e-neighborhood, 
A c (Sf (zo) £ ), carries information on the density of accumulation, as 
well as on the direction of approach to the origin. 

At the end of this Section, we state Proposition 1.3.1. from [S] about 
the formal classification of germs of parabolic diffeomorphisms. We say 
that two germs of diffeomorphisms, f,g& DifftC, 0), are formally con- 
jugate if there exists a formal change of variables <j)(z) = c i zl -> °i ^ 
C, possibly divergent, which transforms / to g, i.e. g(z) = _1 o/o0(z). 

Proposition 1 (Formal classification of parabolic diffeomorphisms, 
Proposition 1.3.1 in [8]). Let 

f(z) = z + ai z k+1 + a 2 z k+2 + ..., ai e C, a* G C, i = 2, 3, . . . , 

fceN, k > 1, 

be a germ of a parabolic diffeomorphism. By a formal change of vari- 
ables, it can be transformed to its formal normal form 

(4) f (z) = z + z k+l + az 2k+1 , aeC. 

The formal type of a diffeomorphism / is thus given by the pair 
(k, a). In the formal change of variables, the multiplicity k remains 
unchanged. The final coefficient a is equal to a = Res( j.A_ z , 0), where 

Res( j( z )_ 2 ? 0) remains unchanged in the formal change of variables, for 
proof see |10| . In literature, this invariant residue is called residual 
index of fixed point zero and denoted by i(f, 0). 

The following remark is important for the sequel. We will see in 
Section [2] that formal changes of variables which are not tangent to the 
identity, i.e. which are of the type (p(z) = C\Z + c 2 z 2 + o(z 3 ), where 
c\ ^ 1, affect the fractal properties, namely the Minkowski content, 
of the diffeomorphism. Therefore, if we want the whole formal class 
of the diffeomorphism /, including its formal normal form fo, to have 
the same fractal properties, we allow only formal changes of variables 
tangent to the identity, i.e. of the type 



00) = z + c 2 z 2 + c 3 z 3 + o(z 3 ). 
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It is easy to check that, if we make formal changes of variables tangent 
to the identity, the first coefficient ai of the diffeomorphism obviously 
remains unchanged, so we have the following version of Proposition [TJ 

Proposition 2. Let f(z), a, k be as in Proposition 1.3.1. By formal 
changes of variables tangent to the identity, f(z) can be transformed to 
its formal normal form 

(5) g (z) = z + a lZ k+1 + a\ ■ a ■ z 2k+1 . 

In this case, the formal type of a diffeomorphism / is given by the 
triple (k, a, aj). 

Note that f is obtained from g if we apply one more change of 
variables <p(z) = Xz, X k = ai, which is not tangent to the identity 
and which eliminates the coefficient a\. To avoid confusion, in the rest 
of the article, the formal normal form Q will be called the standard 
formal normal form and ^ the extended formal normal form, since it 
carries more information on the initial diffeomorphism. 

Let us note that in both cases only first k + 1 coefficients of a diffeo- 
morphism contribute to its formal normal forms. All the monomials 
of order greater than 2k + 1, possibly infinitely many of them, can be 
eliminated one by one by a formal change of variables, without affect- 
ing the former coefficients. Therefore, they are of no importance for 
formal classification, and in the sequel we can restrict ourselves to the 
diffeomorphism up to the order 2k + 1, 

/(*) = z + a lZ k+1 + a 2 z k+2 + ... + a k+l z 2k+l + o(z 2k+1 ). 

Accordingly, to relate formal normal form of a diffeomorphism with 
coefficients in the formal asymptotic development of the directed area 
of the e-neighborhood of an orbit, it will suffice to study only the first 
k + 1 coefficients of the formal asymptotic development. 

2. Main results 

Let us define the directed area and the directed Minkowski content 
of a measurable set. 

Definition 1 (directed area of a measurable set). Let U C C be a 

measurable set, whose center of mass is not the origin. We define the 
directed area of the set U , denoted by A C (U), as the complex number 

A c (U) = A(U)-u t[u) , 

where A(U) denotes the area of U , t(U) G C the center of mass of U 
and Vt{u) — £ C, \ft(u)\ — 1; the normalized center of mass of U . 
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Minkowski content of a measurable set U C C is by definition in Sec- 
tion [T] equal to the the first coefficient in the asymptotic development 
of the area of the e-neighborhood of U. We define directed Minkowski 
content analogously, but using the directed area of the e-neighborhood, 
thus taking into account the position of the set in the plane. 

Let U e denote the e-neighborhood of U. 

Definition 2 (Directed Minkowski content of a measurable set). Let 
U C C be a measurable set, such that its center of mass is not the 
origin. Let d = dim#?7. We define the directed Minkowski content of 
U, denoted by Ai (U), as the complex number 

M C (U) = lim^^. 

Note that, by definition, \M C {U)\ = M(U), where M{U) is the 
Minkowski content of U. 

Let / : C — » C be a germ of a parabolic diffeomorphism: 

f(z) = z + a lZ k+1 + o{z k+1 ), a x e C* , ke N. 

Let the initial point zq lie in an attracting sector of the origin. We 
denote the attracting orbit of / with initial point z by S^(z ). 

In Theorem [3] at the end of this Section, we show that the directed 
area of the e- neighborhood of the orbit S^(zo) has an asymptotic de- 
velopment of the form: 

(6) 

A c (S f (z ) £ ) = K^ 1 ^ + K 2 e l+ ^ + ... + K k ^e 1+ ^ + 

+ K k e 1+ ^ \oge + Se 1+ ^ + K k+1 e 2 \oge + o(e 2 \oge) } e^O. 

Here, coefficients Ki, i = 1, . . . , k + 1, and S are complex numbers. For 
the precise statement and properties of coefficients Ki and S, namely 
their dependence on the coefficients of the diffeomorphism and on the 
initial point, see Theorem [3] at the end of this Section. 

From the development (|6]), it holds that 
A(S f (z ) e ) = \A c (S f (z ) £ )\ = li^+^T +0 (e 1+ wi), e 0. 



Therefore we have that any orbit S^^o) is Minkowski measurable, with: 
(7) 

dim B (^(z )) = 1 - ^-p M(S f (z )) = \K \, M c (S f (z )) = K 1 . 

Motivated by the fact that the first coefficient of ^ incorporates 
directed Minkowski content of the orbit, we define the directed residual 
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content of the orbit as the coefficient in front of the logarithmic term, 
e 2 logs. 

Definition 3 (Directed residual content). We define the directed resid- 
ual content TZ C (S^ (z )) of the orbit S-^(z ) as the complex number 

(8) K c (S f (z )) = K k+1 , 

where Kk+i is the coefficient in front of the logarithmic term e 2 log e in 
the development 

Now we state the two main results of the article. First, the standard 
formal normal form, f(z) = z + z k+l + az 2k+1 , of a given parabolic 
diffeomorphism can be deduced from fractal properties of one of its 
orbits near the origin. 

Theorem 1 (Standard formal normal form and fractal properties of 
an orbit). The standard formal type (k, a) of a parabolic diffeomor- 
phism f(z) is uniquely determined by dimB^S^ " (zq)) , Ai c (S^ (zq)) and 
TZ c (S^(zq)) of any attracting orbit S^(zq) near the origin. 
Moreover, the following explicit formulas hold: 

dim B (S f (z )) 



(9) k 



1 - dim. B (Sf(z )) ' 

fc + 1 _ k+1 . Re ( nC ( Sf ( z °)) \ + 

2 7T V VM c (Sf{z )) J 

+ i-(j)(k)-M(S f (z ))-Im 



n c (Sf(z )) 

l/ M c (Sf{zo)) 

where i / M c {sf(z Q y) is the normalized directed Minkowski content and 4>(k) 
is a function of k, explicitly given by 

_ k(k + 1) J_ F(fT^)+v^ r(l + afei) 

Here, V denotes the gamma function. 

As we have discussed earlier in Section [TJ the converse of Theorem [T] 
is not true. Nevertheless, if we consider the extended formal normal 
form, f(z) = z + a\Z k+1 + a\-a-z 2k+l , instead of standard formal normal 
form, Theorem [T] takes the form of the stronger equivalence statement: 

Theorem 2 (Extended formal normal form and fractal properties of 
an orbit). There exists a bijective correspondence between the following 
triples: 
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(i) the extended formal type of a diffeomorphism, (k, a±, a), 

(ii) (dim B (Sf(z )), M c (Sf(z )), n c (Sf(z ))), 

where S*(zq) is any attracting orbit of a diffeomorphism. The bijective 
correspondence is given by formulas ^ and the following formula for 
a\ : 

(io) a^M%sH^ <- 2 >" < k r 'i + ^-" +1) 



M(Sf(z )) \s/x IXgL) 

The converse states that all the attracting orbits of all the diffeo- 
morphisms of the same extended formal type share the same fractal 
properties. 

Actually, for the precise converse statement, we have to make the 
following remark about the sectorial dependence of fractal properties on 
the initial point of the orbit. Suppose that we know only the extended 
formal type of a diffeomorphism and we want to compute the directed 
Minkowski content and the directed residual content of any attracting 
orbit S^(zq) of the diffeomorphism. The directed Minkowski content is 



given by reformulation of the formula (10): 

where v& is the attracting direction in whose attracting sector zq lies. 
Therefore, the fractal properties do differ slightly in argument for the 
orbits in k different attracting sectors, but they do not differ for the 
orbits inside one attracting sector. Their modules, in particular the 
Minkowski content, are the same in all sectors. 

Remark 1. By ^ and rt8]) ; we can as well express the correspondence 
in Theorem^ and Theorem^in terms of coefficients Ki, K^ + \ and the 
exponent k in the asymptotic development of A (S^(zo) e ), instead in 
terms of fractal properties of the orbit S-'(zo). 



At the end of the Section, we state the auxiliary theorem which 
gives a more precise description of the coefficients in the development 
Q. This theorem is an important part of the proof of Theorem [l] and 
Theorem [2l 

Let f(z) = z + aiz k+1 + a 2 z k+2 + o(z k+2 ), a x 6 O, be a parabolic 
diffeomorphism and let S*(zo) be its orbit with initial point zq lying in 
one of k attracting sectors. Let 



A = (-kai)~* 
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be the one of k attracting directions in whose attracting sector the 
initial condition z$ lies. In other words, we chose the |— th complex 
root of — — whose argument is closest to Zq. By 1/4, we denote the 
normalized complex number A. 

Theorem 3 (Asymptotic development of the directed area of 
^-neighborhoods of orbits). The directed area of e-neighborhood of an 
orbit S-^(zq) has the following asymptotic development: 

(11) 

A c (S f (z )e) = K l£ 1+ ^ + K 2 e 1+ ^ + ... + K k ^e 1+ ^ + 

+ K k e 1+ ^ loge + S{f, z )e 1+ ^r + K k+1 e 2 log e + o{e 2 logs), e 0. 

All coefficients K i} i — 1, . . . , k + 1, are complex-valued functions which 
depend only on k, A and on the first i coefficients a 2 , . . . , of the 
diffeomorphism. The coefficient S(f, zq) is a complex-valued function 
which depends on the whole diffeomorphism f(z) and on the initial 
condition zq. 

Furthermore, 'important' coefficients K\ and K k+ i are of the form: 
(12) 

K k + i r(i + 53 fe) / 2 



Kk+i = va 



k + 1 v a\ 2 ' 

r(|_ 

2fc+2 



2fc+2 I _ rz 

2(fc-l)/|ai|\V(*+i)r(2+ 5 i s ) V7r \ . _ ,a k+1 k + 1 



k + 1 V 2 J r(^) 



+ 



+ g(k,A,a 2 ,...,a k ). 

Here, g(k, A, a 2 , ... , a k ) is a complex-valued function with the property 
g(k,A,0,...,0)=0. 

Note that coefficients K iy i = 1, . . . , k + 1, do not depend on the 
initial point zq, but only on the attracting sector of the initial point 
(via A). Dependence of S on the initial point comes from the directed 
area of the tail of the e-neighborhood of the orbit. On the other hand, 
the first k + 1 coefficients in the development of the directed area of 
the nucleus are independent of the initial point, see Lemmas [2] to [5] in 
Section [3] It is interesting to observe that further coefficients in the 
asymptotic development depend on the initial point zq . 
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3. Proof of Theorem [3] 

We describe below the main steps of the proof of Theorem [3] The 
proof is rather long and technical, so each step is contained in a separate 
lemma below. Some auxiliary propositions are in the Appendix. 

Suppose e > 0. By Definition [TJ 

A c (Sf(z ) £ )=A(Sf(z ) £ )-v t{sf{zo)£) . 

Therefore, we need to compute the first k + 1 terms in the development 
of the area of the ^-neighborhood and the first k + 1 terms in the 
development of its normalized center of mass. Following the idea from 
[13], the e-neighborhood of the orbit, S^(z ) e , can be regarded as a 
disjoint union of the nucleus N £ and the tail T £ . The tail T £ is the union 
of disjoint discs K(zi, e), i — 0, . . . , n e . The nucleus N £ is the union of 
overlapping discs K(zi, e), i — n e +l, . . . , oo. Here, n e denotes the index 
when discs around the points start to overlap. In our case, this 'critical' 
index n £ is unique and well-defined, since the distances between two 
consecutive points are strictly decreasing, see Proposition |2Ji) in the 
Appendix. 

Step 1. In Lemma [TJ we compute the first k + 1 terms in the asymp- 
totic development of the index n £ , as e — > 0. 

Step 2. Using the development for n £ , we compute the first k + 1 
terms in the development of the area of the e-neighborhood of the orbit, 
A(Sf(zo) £ ), as e — > 0. This consists of two parts: first, in Lemma[2j we 
compute the development of the area of the nucleus, A(N E ). Second, in 
Lemma |3j we compute the development of the area of the tail, A(T £ ). 
Finally, 

(13) A(S f (z ) £ ) = A(N £ ) + A(T £ ). 

Step 3. We need to find first k + 1 terms in the development of the 
normalized center of mass of the e-neighborhood of the orbit, ^t(£tf(zo)<0' 
as e 0. Obviously, 

(u) y _ t(N £ )-A(N £ )+t(T £ )-A(T £ ) 

{ ' t(Sf{zoh) |t(iV e ) • A(N £ ) + t(T E ) • A(T £ ) \ ' 

Again, in Lemma |4| we compute first k + 1 terms for the nucleus, 
t(JV e ) • A(N £ ). In Lemma [5j we do the same for the tail, t(T e ) • A(T £ ). 



Now, combining the obtained developments (13) and (14), the result 
for A c (Sf(z ) £ ) follows. 

The following lemmas are used in the proof. They provide asymp- 
totic developments up to the first k + 1 terms of the expressions that 
are neccessary for computing the first k + 1 terms of asymptotic de- 
velopment of the directed area. In all these developments, we provide 
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precise information only on the first and on the (k + l)-st coefficient, 
since they are the only ones that affect the first and (A;+l)-st coefficient 
in the development of the final directed area. 

Lemma 1 (Asymptotic development of n £ ). Suppose n e is the critical 
index separating the nucleus and the tail, as in the proof above. Then 
it has the following asymptotic development: 

(15) 



Pie k + p 2 e 



-i+ 



k+i 



+ p 3 e 1+ s+r + ...+ 



+ p k e + logs + o(loge), e ->■ 0, 



where coefficients Pi = pi(k,A,a 2 , 



k + 1, are real- 



valued functions of k and first i coefficients of f(z). Moreover, 



Pi 



Pk+i 



\\aiA k + l 
-Re 



fc+i 



Qfc+i {k + l)ai 



A 



+ g(k,A;a 2 , . 



.ah) 



ai 2 
. , Ok) is a real-valued function which satisfies 



k + 1 

where g = g(k, A, a 2 
g(k,A,0,...,0)=0. 

Proof By d n = \z n+ i— z n \, n G No, we denote the distances between two 
consecutive points of the orbit. The critical index n £ is then determined 
by the inequalities 

(16) d ns < 2e, d ns ^i > 2e. 

To obtain the asymptotic development of n £ , we first compute as- 
ymptotic development for d n , as n — > oo. Using development (44) for 
z n from Proposition [3] in the Appendix, we get 

(17) 



Zn+l 



aiA 



k+i n -l 



I . _1_2 

k + h 2 n fc +-...+- h 



k™ 



(k + l)a\ 



A 



2k+ 



ik + 1 



+ g(k,A,a 2 , ...,a k ) 



_2_I , 

n k logn+- 



4- o[n h logn), n — > oo, 

where hi = hi(k, A,a 2 , . . . , a^) and g = g(k, A, a 2 , ... , a k ) are complex- 
valued functions and g(k, A, 0, . . . , 0) = 0. Furthermore, 

(18) 



"fc + q 2 n k + . . . + q k n 



k + 1 



aiA k+1 \Re 



a k+1 (k + l)ai 



A*)+r{k,A,a 2 ,..,a k ) 



n 



k log n+ 



on 



-2- 



k logn), n —¥ oo, 
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where = qi(k, A, 0,2, ■ ■ ■ , di) and r = r(k, A, a 2 , . . . , a^) are real-valued 
functions and r(k, A, 0, . . . , 0) = 0. 



From (|16|) and (18) we deduce the asymptotic development of n £ as 



£ — > 0, iteratively, term by term. □ 



Note that the above proof provides developments (17) and (18) for 
z n — z n+\ and for the distances d n between two consecutive points, 
which we also need later. 

Lemma 2 (Asymptotic development of the area of the nucleus). The 

following asymptotic development for the area of the nucleus of the 
e -neighborhood of the orbit holds: 

(19) 

k \Y( 



■ 2 1 2fc+2' 



+ . . . + 4 1} £ 1+ wi + hfe 2 log £ + o(e 2 loge), e -> 0. 

Here, hi — hi(k, A, . . . , Oj) , i = 2, . . . , k, are real-valued functions of 
k and first i coefficients of f(z). 



Proof. By Proposition |4|n) in the Appendix, the area of the nucleus 
can be computed by adding areas of infinitely many crescent-shaped 
contributions. Furthermore, Proposition [5] provides the formula for 
computing such areas. We have 




(20) A(N £ ) = 2e 2 V I = d 1 - || + arcsin dr ' 



x . As 2 2e 

n=n e +l \ , 

By Proposition [6] in the Appendix, this sum can be replaced by the 

following integral: 

(21) 

A{N £ ) = 2e 2 I + arcsin^ | dx+Otf), £ -+ 0, 



where d(x) is the strictly decreasing function from Proposition |6j 

d(x) = qix^ 1 ^^ -f g 2 x~ 1_ fc + . . . + qkX~ 2 + qk+\X~ 2 ~k logs + Dx~ 2 ~k . 

We now compute the first fc + 1 terms in the asymptotic development 



of the integral from (21), as e — > 0. Applying the change of variables 
t = we get 



(22) I = -2e 2 " (tVl - 1 2 + arcsin t) 1 dt. 

Jq V / d'{x(t)) 
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Here, x[t) = d~ l {2et). Note that, for a given e, t is bounded in [0, 1). 
Therefore it holds that: 

(23) (et) — > 0, as £ — > 0, uniformly in t. 

The development of x(t) = c/ _1 (2et), as e — > 0, can be deduced using 
the already computed development for n £ = d~ 1 (2e) in Lemma [TJ We 
have that 

(24) 

1 k 



-2- 2+ ^T|a 1 A fc+1 | 1 ^(^)- 2+ ^i+p 2 (et)- 2+ ^T + 
d'(x(t)) k + 1 

+ . . .+p{ 1) (ety 2+ ^T +pf(et)- 1 log(et) + 0((rf)- 1+ ^i), 



2s 



where pi = pi(k, A,a2, ■ ■ ■ ,ai), i = 2, . . . , k, are real- valued functions. 

Using Proposition [7|in the Appendix, we remove £ from the boundary 
of /. The integral in (|22j) is equal to 



r 1 i 

(25) I = -2e (Wl-t 2 + arcsin t) , . . . . dt + o(log < 

y d'(x(t)) 



s 0. 



Now, substituting the development (24) in (|25j) and using (23) to 
evaluate the last term, we get 

(26) 

( 2 Y 2+ ^ k _ , . i m , , 2 

'=(]MFi|J jfTT -T 1 . e -'+™-2 ft .r,. e -^- 

- ... - 2p<" ■ T t ■ _ 2 p' 2) ■ S M • loge + o(loge), e -> 0. 



Here, functions are real- valued functions from the development (24) 
and Sf. + i and T», i = 1, . . . , k, are the following finite integrals: 

Ti= f (ty/l - t 2 + arcsin t)t~ 2+ ^dt, i — l,...,k, 
Jo 

S k+1 = (ty/l -t 2 + arcsin t)t _1 logtdfc 

Since Ti = f JilSi, - A) , combining (pB, and ph, 

V 1 \ 2 + 2fc+ 2 J / 

we get the development (19) for A(N £ ). □ 

Lemma 3 (Asymptotic development of the area of the tail). The area 
of the tail of the e -neighborhood of the orbit has the following asymptotic 
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A(T £ ) =n(^—Y 1+kll z 1+ ^ + f 2 s 1+ ^ + ... + f k s 1+ ^ + 



e 2 log e+ 



development: 
(27) 

{\a 1 A k + 1 r 

+ o(e 2 logs), e -»• 0. 

ifere, /j = A, a 2 , . . . , a*), z = 2, . . . , k, are real-valued functions 
which depend only on k and the first i coefficients of f(z) . The function 
g = g(k, A, a 2 , . . . , Ojj.) has the property that g(k, A, 0, . . . , 0) = 0. 

Proof. Since the tail, by definition, consists of n £ — 1 disjoint e-discs, 
we have that \T £ \ = (n e — 1) ■ e 2 it. The statement follows from ( 15 ). □ 



Lemma 4 (Asymptotic development of the center of mass of the nu- 
cleus). Let t(N e ) denote the center of mass of the nucleus of the e- 
neighborhood. The following asymptotic development holds: 

(28) t{N £ ) ■ A(N E ) = qi e 1+ ^ + q 2 e 1+ ^ + q 3 e 1+ ^ + ... + q k e 2 + 

+ qk+i£ 2+1 ^ loge + o(e 2+£ TT loge), e ->■ 0. 

Here, qi = qi(k,A,a 2 ,...,cii), i = 1, . . . , k + 1, are complex-valued 
functions which depend on k and on the first i coefficients of f(z). 
More precisely, 



u 2 



91 2^-1)^(1 + ^) ^JIk^iJ ' a 
_ / r _ ni+jM. \ ( 2 



i 

fe+i 



Jm (-¥ - —J— ) -A-i + g(k, A, a 2 , . . . , a fe ), 



where g = g(k, A,a 2 , ■ ■ ■ , afc) a complex-valued function such that 
g(k,A,0,...,0)=0. 
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Proof. By definition of the centre of mass and by Propositions |4ju) and 
[5] in the Appendix, we have that 



n=n £ +l 



t(iV £ ) • A(N £ ) = z nE ■ e 2 n + A(D n )t(D n 

00 

+ E 



dm 



n=n e +l 



2s 



d 2 . 



4e 2 



+ arcsm 



2e 



Zn+ 



+- E 

n=n e +l 



r/ 2 , 



, , 1 — -^r — arcsin W 1 
2e V 4e 2 V 4e 2 



r/ 2 



J ( z n ~ z n+l)- 



Here, D n , n > n e , denote the contributions to the nucleus from the 
£-discs of the points z n . 

We first show that 
(29) 



t(JV e ) ■ A(N, 



E 



n=n e + l 



+ arcsin 



z n + 0(e 2+ ^i' 



2e V " 4e 2 ' 2e 
as e -» 0. From ((To} and rt44}, 2 ne ■ e 2 n = 0{e 2+ ^+i), as e -)> 0. 



On the other hand, by (17), we have that z n — z n+ i = C^n -- * 1 ), as 
n — > 00. Therefore, using boundedness of the term in parenthesis, 
integral approximation of the sum and then (15), we get 



E 

n=n £ +l 



dn / 1 d n . / d„ 

— \ 1 — ~i — o — arcsm \ 1 - 
2e\ As 2 V 4e 



<^ 2 E 



^"J (^n — ^n+l) < 



n=n e +l 



for some constant C > 0. This proves (29). 



To compute the first k + 1 terms in the asymptotic development of 

d n 



the sum in (29) 
(30) 



S 



E 

n=n e +l 



2s 



1 



4s 2 



+ arcsin 



dn 
2s 



as s — » 0, we use the same idea as in Lemma [2] Therefore we omit 
the details. To make the integral approximation of the sum S, we have 
to cut off the formal developments d n and z n to finitely many terms. 
Let d* n be as in Proposition |6j d n = Jk+id n + Dn~ 2 ~k. By Jk+iz n , we 
denote the first k + 1 terms in the asymptotic development of z n . It 
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can be shown similarly as before that 



S 



n=n £ +l 



2e 



^g- + arcsm ^ ) J k z r , 



+ o(e k + 1 \oge). 



Since the real and the imaginary part of the function under the summa- 
tion sign are strictly decreasing, as n — > oo, we can make the integral 
approximation of the sum: 



(31) S 



d(x) 
~2e 



d(x) 2 . d(x) 

+ arcsm 



Ae 2 



2e 



z(x)dx + o(e k + 1 loge). 



The function d(x) is as defined in (54) of the Appendix, and z(x) is 
equal to 
(32) 

z(x) = g x x~ k +g 2 x~ k + g 3 x' k + 9i x' k + . . . + g k x 



-i 



■g k+ ix * logx 



with coefficients gi e C from the development (44) of z n in the Appen- 
dix. 

By making the change of variables t = 4^ in the integral, we get 
(33) 

S = -2e 
as e — > 0. 



l+0(e S+T) 



(Wl - t 2 + arcsint^ Z jf^). dt + o(£ fe +i loge) 
V / d'(x(t)) 



Using (24), (32) and the development for x(t) from the proof of 
Lemma 21 after some computation we get the development for , 
as et — >■ T). Again, let us note that e£ — >■ uniformly in t, as £ — > 0, see 
(23) before. Substituting the development in (33) and proceeding in a 
similar way as in Lemma |5J we get the development ( 28 ) . □ 

Lemma 5 (Development of the center of mass of the tail). The fol- 
lowing development for the center of the mass of the tail of the e- 
neighborhood holds: 

(34) 

t{T £ )A{T £ ) 



k-1 



-7T 



7T 



fc + 1 V |oi^ fc+1 | 



i 

fc+i 



k + l 



■ A ■ e 1+ fc+i + g 2 e 1+ + 

+ ...+ 5fc_ie 1+ sri + g k e 2 log e + 5(/, z )e 2 + 

Im(^±I - ■ i ■ A + A, a 2 , . . . , a fc ) 

aj 2 

• e 2+ fcTT loge + o(e 2+ fcTT loge), e -> 0. 
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Here, gi = gi(k,A,a,2, ■ ■ ■ ,a>i), i = 2,...,k, are complex-valued func- 
tions of k, A and first i coefficients of f(z). The function S(f,zo) is 
a complex-valued function which depends on the whole f(z) and on the 
initial condition zq. The function h = h(k, A, 02, ... , a&) is a complex- 
valued function which satisfies h(k, A, 0, . . . , 0) =0. 

Proof. 

t(T £ ) • A(T £ ) = ^ n = l A( Z T n ; £n ■ A(T £ ) = e\ £ = 

A ^ e > n=l 

n E — l 

= E 2 7l(z + ... 2„(/, 20 )) + S 2 7l ■ ^ Z n- 

n=n(f,z ) 

Here, n(f, zq) is chosen to be the first index, obviously depending on 
the diffeomorphism / and on the initial condition zq, such that 

z n = Jk+\z n + R(n), where < Cn _1 ~£, for n > n(f, z ), 

for some constant C > 0. Then 
(35) 

n E — 1 n e — 1 n E - 1 n E — l 

^2 z n = 9i ^2 n ~ % +92 n~* + ...+g k ^ n ~ 1+ 

n=n(f,z ) n=n(f,z ) n=n(f,z ) n=n(f,z ) 

n e — 1 n e — 1 

+ 9k+i n~ 1 ~*logn + 5^ ^( n )' 

n=n(f,z ) n=n(f,z ) 

where complex numbers ^ are as in the development of z n , see Propo- 
sition [3] in the Appendix. 

We now compute the first k+1 terms in the asymptotic developments 
of (35), as n £ — » oo. 



Firstly, we concentrate on the last sum in (35). We show that 

n 

(36) Yl R ^ = C ^ /) + °( n ~*)> " °°» 

l=n(z ,f) 

where #(/) = 0(r 1_ «), as / — > oo, and C(z , /) is a complex constant 
depending on the diffeomorphism and on the initial condition. From 
the asymptotics of R(l), the sum Ylu= n (z /) -^(0 * s obviously convergent 
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and equal to some constant C(zq, f) G IR. We write 

n oo oo 

E m)= E a(0-X>w = 

l=n(zo,f) l=n(zo,f) l=n 

= C(z , f) + 0(n~«), n -> oo, 

where the second sum is evaluated as 0(n _1 ' fe ) by integral approxima- 
tion of the sum. 

Secondly, we estimate first three terms in the asymptotic develop- 
ments of the first k + 1 sums in (35), as ne — > oo. We show the 
procedure on the first sum. Let 

n 

F(n)= E 

l=n(f,z ) 

Obviously, it satisfies the recurrence relation 

F(n + 1) - F(n) = (n + l) - *, nGN, 

with initial condition F(n(f, Zq)) — n(f, Zo)~k . We determine the first 
term in its development by integral approximation: 

F(n) = T^n^r + 

/b 1 



fe-1 . 



where i?(n) = o(n~^) as n — > oo. Using this development, from 
recurrence relation for F(n) we get the recurrence relation for R(n) 
and the initial condition R(n(zo, /)). By recursion, we get 

n 

R(n) = R(n(z , /)) + E ^n 1 ^). 

'="(20,/) 

Using ((36]), we conclude that Yl?=n(f,z ) ^ = C{n(z ,f)) +0(n~ 1/k ), 
n — > oo. The same procedure can be repeated for other sums. 

Thus we obtain the development of the sum Y^= n (zo /) Zn -> as n£ ~~ 
oo. Substituting n e with the development (15), we get the development 
(g, as e ->• 0. □ 

Proof of Theorem [3| The proof follows from Lemmas [T] to |5j as 
described at the beginning of the Section. □ 

4. Proof of Theorems Q] and [2] 

We now prove the main results, Theorem [T] and Theorem |2j In the 
proofs, we need the following lemma. It shows that the leading expo- 
nent and the relevant first and (k + l)-st coefficient in the development 
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of the directed area remain unchanged by a change of variables tangent 
to the identity, transforming the diffeomorphism to its extended formal 
normal form. 

Lemma 6 (Invariance of fractal properties in the extended formal 
class). Let f\(z) and /^(-s) be two germs of parabolic diffeomorphisms 
which belong to the same extended formal class (k, a, a{). Then it holds: 

&m B {S h {w Q ))=d\m B {S h {z Q )), 

M c (S fl (w ))=M c (S f *(z )), 

n c (s^(w )) = n c (s h (z )). 

Here, zq and wo are any two initial points chosen from the attracting 
sectors of f\ and f2 with the same attracting direction. 

In the proof of Lemma [6] we need the following auxiliary lemma: 

Lemma 7. Let f(z) be a parabolic diffeomorphism and let g(z) = 
o / o <j>i(z), where 4>i(z) = z + cz l , I > 2. Let S^(z ) = {z n } 
be an attracting orbit of f(z) and let S 9 (v ) = {w n = (pi(z n )} be the 
corresponding attracting orbit of g(z). Then it holds that 

(37) K? {zo) = K 8 ^, = K S k ^°\ 

where K\ and K k+ i denote the first and the (k + l)-st coefficients 



in the asymptotic developments (11) of the directed areas of the e 



neighborhoods of the corresponding orbits. Furthermore, the equalities 



(37) hold also if S^(z ) and S 9 (vq) are any two orbits of f(z) and g{z) 



respectively which converge to the same attracting direction. 

Proof. Let {z n } be an attracting orbit of f(z). We first take S 9 (vq) = 
{w n } to be the image of {z n } under (pi, I > 1. Using development 



(44) for z n , we compute the development of w n = 4>i(z n ). It is easy to 
see that, since / > 1, the first coefficient and the (k + l)-st coefficient 
remain the same as in z n , while the other coefficients can change. In 
particular, the attracting direction A for S 9 (vq) remains the same as 
forS'(zo). 

On the other hand, it can be seen in Section [3] that only the first 
and the (k + l)-st coefficient of the development of z n participate in 
the first and the (k + l)-st coefficient of the developments of z n — 
z n+ \, d n and n £ . Finally, the first and the (k + l)-st coefficient in the 
development of A C (S^ (zo) e ), K\ and Kk+i, depend only on the first 
and the (k + l)-st coefficient in the development of z n , not on other 
coefficients. Therefore, the two coefficients remain unchanged in the 
change of variables 4>i(z). Finally, since K 1 and K k+ i do not depend 



e-NEIGHBORHOODS OF ORBITS OF PARABOLIC DIFFEOMORPHISMS 21 



on the choice of initial point zq and vq inside one sector, we can choose 
any two orbits of the initial and of the transformed diffeomorphism 
converging to the same attracting direction A. □ 

Proof of LemmauA For the given diffeomorphisms fi(z) and f2(z), let 

1 o 12 12 j 2 

4> = (fik ^fe'-i ° ■ ■ ■°4>2 denote the changes of variables obtained by 
composition of k — 1 transformations of the above type, which present 
the first k — 1 steps in transforming fx and f2 to their extended formal 
normal forms. Let 

(38) g x =(0 1 )" 1 of 1 o<f) 1 =z + a x z k+l + a\az 2k+l + ..., 
92 =(^-i o / 2 o 2 = * + + a 2 az 2fc+1 + . . . . 

Obviously, by Lemma [7| it holds that: 

(39) Kf = K{\ Kf = Kf and = K* v K» t = 

for the orbits corresponding to the same attracting direction. The 
notation K( is a bit imprecise, since the value differs for orbits in k 
sectors, but we use it for simplicity and keep in mind that we always 
consider orbits converging to the same attracting direction. 
Let go be the extended formal normal form, go(z) = z + aiz k+1 + 
a\az 2k+1 . By further changes of variables, transforming gi and g2 to 



the extended formal normal form (fa; the (2/c + l)-jets from (38) remain 
the same. Therefore we have, by the development (12) in Theorem [3j 
that 

(40) iff = iff, Kf = Kf and K°} +1 = Kf +V K» x = K 9 k ° +1 , 



for the orbits corresponding to the same attracting direction. By (39) 
and gjf, it follows that K{ x = K( 2 and = K f k 2 +1 , for the orbits 



of fx and f2 converging to the same attracting direction. 

Finally, changes of variables do not change the multiplicity + 1 
of the diffeomorphism. Therefore the leading exponent of the directed 
areas for all the orbits equals 1 — k ~^. 

Relating the coefficients K%, Kk+i and exponent k with fractal prop- 
erties of orbits, by ^ and ([8]), the statement follows. □ 

Note that the statement of the above Lemma is no longer true if we 
admit changes of variables which are not tangent to the identity. Only 
box dimension is then preserved. 



Proof of Theorem^ Let f(z) = z + a\z k+1 + o(z k+1 ) be a parabolic 
germ and let go(z) = z + aiz k+1 + a\ ■ a ■ z 2k+1 be its extended formal 
normal form. Let S^(z ) t> e an attracting orbit of f(z) and let S 90 (w ) 
be an attracting orbit of go(z) with the same attracting direction. 
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The bijective correspondence between k and dini^S'-' (zo)) is obvious 
by ([7]). Let k then be fixed. Applying formulas (12) from Theorem [3] to 
the orbit of the formal normal form go(z), we get the following formulas: 

(41) 



K{ 



K 90 



fc + 1 



2(k-l) 
fc + 1 



r(i 



1 



2k+2 > 



1 



2A:+2< 



7T 



N 
2 



T ite(< 



2 

! a il 
fc + i 



l/(*+l) 



) + 



'2 T 2fc+2 j 

i/(k+i) r(2 +55 ^) 



— \/7T 



• 1 



Im(a — ) 



K 90 



By Lemma [6j 

(42) K{ = K 9 \ K[ +1 
On the other hand, by |7| and ([8]), 

(43) tff = -M c (^(zo)), < x = n c (Sf(z )). 

Using (42) and (43), we see that formulas ^ and (10) in Theorem|2]are 



just reformulations of (41). They give, for a fixed k, the bijective cor- 



respondence between the pairs (ai,a) and (Ai c (S^ (z )), TZ c (S^(z ))). 

□ 

Proof of Theorem^ Let f(z) and go{z) be as in the above proof. The 
standard formal normal form fo{z) is given by fo(z) = z + z k + az 2k+l , 
where a is the same as in the extended form go{z). The normal form 
fo(z) is obtained from g (z) by making one extra change of variables 
of the type 

cj)(z) = a l 1 ' k z, 

in order to make coefficient a\ equal to 1. Since a and k are the same 
as in go(z), formulas ^ expressing k and a from the fractal prop- 
erties of the orbit S* (zq) have already been obtained in the proof of 
Theorem [U Therefore the standard formal normal form of a diffeo- 
morphism, described by the pair (k,a), can be deduced from fractal 
properties {dim B (S f {z )), M c {S f {z ), K c (S f (z ))) of just one orbit 
of the diffeomorphism. □ 

Let us note that Theorem [T] cannot be formulated as an equivalence 
statement between (k,a) and fractal properties. From the pair (k,a), 
one cannot uniquely determine the fractal properties of the orbit of 
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the initial diffeomorphism f(z). Aside from the box dimension, the 
diffeomorphisms from the same standard formal class do not share the 
same fractal properties. By ^ in Theorem |3j the directed Minkowski 
and residual content depend on the first coefficient a±. The information 
on the initial fractal properties is lost by making changes of variables 
which are not tangent to the identity and which change a±. 

5. Perspectives 
Let us comment shortly on the perspectives for further research. 

5.1. Problem of analytic classification. 

We have shown in this article that the formal type of a diffeomor- 
phism can be read from any orbit, using only its fractal properties. We 
are interested in prospects of analytic classification of parabolic diffeo- 
morphisms using e -neighborhoods of orbits. The analytic classification 
was given independently by Ecalle and Voronin in [2] and [15]. The 
analytic classes are given by the formal invariants (k, a), as well as by 
2k diffeomorphisms, which are called Ecalle- Voronin functional moduli 
of analytic classification. 

In general, a diffeomorphism is analitically conjugate to its formal 
normal form only sectorially. One orbit of a diffeomorphism lies com- 
pletely in one sector. Our goal is to see if the analytic type can be 
read from the directed area of the e-neighborhoods of only one orbit, 
or if perhaps something can be said in this directon if we consider e- 
neigborhoods of one orbit per sector for 2k sectors, and compare them 
in an appropriate way. 

5.2. Can one hear the shape of a drum? In this article we were 
motivated by the question: to what extent a parabolic diffeomorphism 
itself can be reconstructed from its one realization, that is, from its one 
orbit? So far, we know that, from only one orbit, we can tell the formal 
type of a diffeomorphism. The concepts are somewhat similar to the 
concepts of the famous problem: Can one hear a shape of a drum?, 
presented by M. Kac in 1966. The question that is posed is if one can 
reconstruct the equation from only one solution, or, if not completely, 
how much can be said. 

The vibrations of a drum are given by the Laplace equation with zero 
boundary condition on a given domain Q. The domain of the equation 
is the only unknown in the problem. The eigenvalues of the Laplace 
operator, < Ai < A 2 < • • ., A, — > oo, present the frequencies. They 
are coefficients in the Fourier development of the solution. One tries 
to reconstruct the domain of the equation from these eigenvalues. 
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Let N(X) = {\i : Aj < A} be the eigenvalue counting function for the 
Laplace operator on Q. It was conjectured that from the asymptotic 
development of N(X), as A — > oo, one can obtain some properties of 
the domain: 

Conjecture (Modified Weyl-Berry conjecture, Conjecture 5.1 in [7j). 
Iffl C R has a Minkowski measurable boundary T, with box dimension 
de(N - 1,N), then 

N(X) = {2uY N B N ■ A(Q)X N/2 + c Ntd ■ MiT) ■ X% + o(A?), A ^ oo. 

Here, Bn is the volume of the unit ball in WL N , A(Q) the Lebesgue 
measure of the set Q G M. N and -M(T) the Minkowski content of the 
boundary. The constant CN,d is a real constant depending only on N 
and d. 

The conjecture was proven in the one- dimensional case, N = 1, in 
Corollary 2.3 in [7]. In other dimensions, it is still open. 

Although we do not see the direct relation between two problems, in 
many aspects they appear similar. The general idea of reconstructing 
the equation from one solution is common, as well as the fact that, for 
obtaining more information on the equation, we need to utilize further 
terms in appropriate developments. 

For more on this problem, see [B] or [7]. 



6. Appendix 

Here we state auxiliary propositions that we need in the proof of 
Theorem HI 

Let f(z) = z + ai z k+1 + a 2 z k+2 + a 3 z k+3 + . . ., a { G C, a 1 ^ 0, 
be a parabolic diffeomorphism. Let the initial point Zq belong to an 
attracting sector. We denote by A the attracting direction 

A = (-kai)~^, 

where we chose the one of k complex roots for which z is closest to 
the direction A. 

Proposition 3 (Asymptotic development of z n ). Let z n = f( on ^(zo), 
n G No, denote the points of the orbit S^(zo)- Then 

(44) 



gin k + g 2 n k + g 3 n >° + g 4 n *+... + gun + 

+ gu+in k logn + o(n k iognj 



e-NEIGHBORHOODS OF ORBITS OF PARABOLIC DIFFEOMORPHISMS 25 

where coefficients gi = gi(k, A, ci2, ... , ai), i = 1, . . . , k + 1, are complex- 
valued functions of k and first i coefficients of f(z). More precisely, 

9i = A, g k+l = -\a^ - + h(k, A,a 2 ,..., a,)) , 

where h = h(k, A, a 2 , ... , a k ) is a complex-valued function which satis- 
fies h(k, A, 0, . . . , 0) =0. 

Proof. The following proof mimics the technique for obtaining the as- 
ymptotic development of a real iterative sequence from [TJ, Chapter 
8.4. In the complex case, we repeat the whole technique sectorially. 
Suppose as above that z lies in an attracting sector around attracting 
vector A. By [TO], the whole trajectory {z n } lies in that attracting 
sector and is tangent to A at the origin. On this sector, the change of 
variables 



i 



(45) z = Aw fc 

is well-defined, the complex root of w being uniquely determined. The 
trajectory {z n } is transformed to {w n } and obviously 

_ i 

(46) Arg(w n k ) — > 0, as n — > oo. 
The recurrence relation for z n 

(47) z n+1 = z n + a lZ k n +l + a 2 z k n +2 + a 3 z k+3 + ... 
transforms to the following recurrence relation for w n : 
(48) 



A k w~ l + o( W - lN 

ai V «i 2 J 

Obviously, = iY.tM ~ w i-i)- B y @> ^ holds that ( w i ~ 

Wi-i) — y 1, as I — > oo, therefore lim n ^oo — = 1. From (48) we then 
have 

w n+ i -w n = l + 0{n~k). 
By recursion and using integral approximation of the sum, we get 

w n = n + 0{n~). 
To compute the exact constant of the second term, with this develop- 



W n+ i =w n 


+ 1 + 


a 2 -i 
ai 


ai 


+ 


^A k ~ 


fe-i 

* + 





ment, we return to (48) and get 



2 _I „/ _2 

w n+ i - u; n = 1 H An * + Oin k 

a 
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By recursion and using integral approximation of the sum, 

. k k=± , k-2 

w n — n-\ A- n k 4- (j[n * j. 

a k — 1 

Repeating this procedure k times, we get the first k + 1 terms in the 
development of w n as n — » oo. The development for z n then follows 



from the development for w n , (45) and (46). □ 



The next two propositions give the tool for computing areas and cen- 
ters of mass of e- neighborhoods of orbits. Let f(z) be a parabolic diffeo- 
morphism and S^(zq) = {z n , n G No} its attracting orbit. Let K(zi,e) 
denote the e-disc centered at zi. We represent the e-neighborhood of 

Sf(z ) as 

oo 

(49) S f (z ) £ = \jD l . 

i=0 

Here, D = K(z ,e) and Di = K(zi,£)\\J^l K(zj,e), i £ N, are 
contributions from £-discs of points z^. 

Proposition 4 (Geometry of ^-neighborhoods of orbits). 

(i) Distances between two consecutive points of the orbit, \z n+ i—z n \, 
are, starting from some n^, strictly decreasing as n — > oo . 

(ii) For small enough e > 0, 

i-l 

K( Zi , e)\ |J K( Zj , e) = K{z h e)\K{ Zi ^e), i e N. 

Proposition |3J(m) means that all contributions Di are in crescent or 
full-disc form, determined only by the distance to the previous point 
Zi-i. The positions of the points zq, . . . , Zi-2 do not affect the shape of 
Di, see Figure |6l 



Proof, (i) Using development (17), we compute that 



k + 1 2fc+l 2fc+l 

z n - z n+l = (z n+1 - z n+2 ) + A—^-n k + {n * ) 

and 

A _fc+i fc+i 
z n+1 - z n+2 = j-n fc +o(n «= ). 

Obviously, in the limit as n — )■ oo, their argument is the same. For n 
big enough, the value of the nonordered angle between z n — z n+ i and 
z n+ \ — z n+ 2 is therefore less than |, so \z n — z n+ i\ > \z n+ i — z n+ 2\ for 
n big enough. 




FIGURE 2. 1. admissible position of discs, 2. nonadmissible 
position of discs. 

(ii) Let T n denote the midpoint and s n the bisector of the segment 
[2 n+ i,z n ], n G N. It will suffice to show that there exists e > such 
that for every n > no and for every k G N, the distance from the 
intersection of s n and s n+ k, denoted S nt k, to the midpoint T n is greater 
than e. In this way we ensure that the union of intersections of e-disc 
of each new point of the orbit with the £-discs of all the previous points 
is a subset of the intersection with the £-disc of the previous point only. 

We first show that the two consecutive bisectors s n and s n+ i intersect 
at the distance from T n which is bounded from below by a positive 
constant, as n — > oo. 

The bisector s n can obviously be parametrized as follows 

(50) T n + t ■ i(zn - z n+1 ) = Zn + ^ n+1 + t ■ i(z n - z n+ i), t G R. 

We denote by t n G R the parameter of the intersection S n \ of s n and 
s n+ i. The complex number 

z n T . . •/ \ rp Z n z n+2 , •/ \ 

2 r* t-n ' 1\Zn z n+l) -'n+l ^ n ' H^i z n+l) 

is perpendicular to ^ n+1 — 2;„ +2 . Therefore their scalar product, denoted 
by (.|.), is equal to 0, and we get: 

, _1 (z n — Z n+ i\z n+ i — Zn+2) + \z n +l — Z n +2\ 2 

2 Re(z n - z n+1 )Im(z n+ i - z n+2 ) - Im(z n - z n+1 )Re(z n+1 - z n+2 ) ' 



Using development (17), after some computation, we get that the de- 

• /^/ _3fc£3T . 3|AI 2 _2fc+2 . 2fc+2 . 

nominator is 0(n * ), while the numerator is js-n k +o{n fc ). 

fc+ 1 

Therefore, t n > Cn~k~ , for some positive constant C > and n > n . 
Since \z n — z n+ i\ ~ n~~ , the distance 

d(T n , S n ^i) |t n | • \z n ^n+l I 

is bounded from below by some positive constant for n > n , say by 
M > 0. 
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It is left to show that the same lower bound holds not only for con- 
secutive, but for any two bisectors s n and s n+ k, k G N, n > n . We can 
see from the development (17) that the points of the orbit approach 
the origin in the direction A. We draw the stripe of width M/2 on 
both sides of that tangent direction. Obviously, for n big enough, no 
two bisectors can intersect inside the stripe without two consecutive 
bisectors being intersected inside the stripe, which is a contradiction 
with the first part. Therefore, the distances from the midpoints to the 
intersections of corresponding bisectors when n — > oo are uniformly 
bounded from below by e.g. M/4. 

Taking e < M/4, we have proven the statement. □ 

Proposition 5. Let z, w G C (or IR 2 ), e > 0. Suppose \z — w\ < 2e. 
Let D denote the crescent D = K(z,e)\K(w,e). Then its area is equal 
to 



\z — w\ L \z — w\ 2 . \z — w\ 



(01, ^(D) = 2:-(^^yi-^-^ + arcsin } 
and its center of mass is equal to 



1^1 J I - - arcsin J I - 
(52) t(D) = z + e 2 (w-z)^ V 

Proof. The proposition is proved by integration, 

M D ) = J J dxd V, tp) = ~^Jyj j xdxdy +i- J J ydxdy 

□ 

The following two propositions are auxiliary results in the proof of 
Lemma |2j 

Proposition 6. The sum 



is equal to 



d(x) I d(x) 2 . d(x) . , /A , , 

2TV 4^ + arcsm 2T ,( + ( K 

as e — > 0, where d(x) is given by 
(54) 

d(x) = qix~ x ~* + q 2 x~ 1 ~h + . . . + q k x~ 2 + q k +ix~ 2 ~k \ Q g X + Dx~ 2 
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All the coefficients qi are the same as in development (18) of d n and 
D Gl is some constant. 

Proof. The idea is to apply integral approximation of the sum. The 
problem is that we only have formal asymptotic development of d n . 
The idea is to cut off the formal asymptotic development at (k + 1)- 
st term, to get a continuous and decreasing function of n under the 
summation sign. We show here that the cut-off remainder is in some 
sense small and contributes to the sum with no more than 0(1), as 
e -> 0. 

We denote by Jk+id n the first k + 1 terms in the asymptotic expan- 



sion (18). For the sum with truncated d n , 



oo 



(55) I Jj ^r V 1 _ {Jk ll d 2 n)2 + arcsin Jj ^ n 



n=n E 



to be well-defined, we have to ensure that < Jk+\d n < 2e for n > 



n e . Since d n < 2e for n > n £ by (16), it is enough to achieve that 
Jk+\d n < d n , for n > n £ , where e is sufficiently small. This is obtained 
by adding the term Dn~ 2 ~k to Jk+\d n . Here, D is chosen negative and 
sufficiently big by absolute value. We denote d* n = Jk+id n + Dn~ 2 ~k. 
Obviously, 

(56) d n = d* n + 0(n- 2 - 1 *). 



Let us denote the function under the summation sign in (53) by h(x): 



XI X 2 X 



h(x) = -^l- — + arcsin ^ 



Then, h'(x) = 1 — (^) . By (56) and by the mean value theorem, 

(57) h(d n ) = h(d* n ) + h'(£ n ) ■ 0(n- 2 -^), e„ e K, d n \. 
Furthermore, 

(58) < h'{i n ) <-, n> n £ . 



The initial sum (|53j) can, by (|57j), be evaluated as follows: 

oo oo oo 

(59) S = ]T h(d n ) = Kd*n) + J2 h'(QO(n- 2 -- k ). 

n=n E n=n E n=n E 

By ( J58| and Lemma [TJ using integral approximation of the sum, we get 

oo f oo f 

(60) I ]T h'(U ■ 0{n- 2 -l)\ <yY, n ~ 2 ~ l < < C ' 
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for some constant C > 0, as e — > 0. 

Furthermore, using integral approximation of the sum and the fact that 

the subintegral function is bounded from above, we get 

(61) 

d(x 



2s 



1 



d(x) 2 d(x) 
+ arcsin 



As 2 



2s 



Finally, by (59), (60) and (61), the result follows. 
Proposition 7. The integral 

. d(n e ) 



dx+0(l), e^O. 

□ 



/ " ( Wi - t 2 + arcsin t ) — — -—-dt, 
Jo V > d\x(t)) 

where d(x), x(t) and n £ are as in Lemma^ is equal to the integral 

/ (ty/l — t 2 + arcsin t) -, dt + o(£ _1 ), 
J \ / d'[x(t)) 



as e — > 0. 



Proof. We first show that the upper boundary in the integral is 
equal to 

/ >. d(n £ 
(62) ±± 



1 + 0(e 1 ~^r), e^0. 



By (jig and (g§, 
(63) 



d(n £ ) = d* ne = d ns + 0{e" 



fc+i 



From (18), it can easily be seen that d n+ \ = d n + 0(n fe), thus by 
(15) and (16), we get 



(64) 



d nE = 2e + 0(e 



2- 



fe+i 



)• 



Combining (63) and (64), (62) follows. 



. d(n e ) 



Using (62), the above integral / = f Q 2e (tyl — t 2 + arcsint) d ,^ t ^ dt 
can be written as the sum 



/■i i 

I— (ty 1 — t 2 + arcsin t) - ,, dt+ 
Jo d'{x{t)) 



+ 



l+0(e fc +!) 



(tyl — t 2 + arcsin t) — 



d'(x(t)) 



dt. 



By (24), 



i 



d'(x(t)) 



0((et) 2+k + 1 ). It is then easy to see that the sec- 



ond integral equals 0{e : ), as e — > 0, due to the boundedness of the 
subintegral function in the neighborhood of t = 1. □ 
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